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We consider quantum mechanics on constrained surfaces which have non-Euclidean metrics and 
variable Gaussian curvature. The old controversy about the ambiguities involving terms in the 
QO [ Hamiltonian of order Ti 2 multiplying the Gaussian curvature is addressed. We set out to clarify 

■ the matter by considering constraints to be the limits of large restoring forces as the constraint 

coordinates deviate from their constrained values. We find additional ambiguous terms of order 
h 2 involving freedom in the constraining potentials, demonstrating that the classical constrained 
Hamiltonian or Lagrangian cannot uniquely specify the quantization: the ambiguity of directly 
quantizing a constrained system is inherently unresolvable. However, there is never any problem with 
a physical quantum system, which cannot have infinite constraint forces and always fluctuates around 
the mean constraint values. The issue is addressed from the perspectives of adiabatic approximations 
in quantum mechanics, Feynman path integrals, and semiclassically in terms of adiabatic actions. 
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> ! I. INTRODUCTION 

l> 
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The controversy over the proper way to quantize a constrained dynamical system has had a long and interesting 
history in twentieth century physics. One treatment of the subject opens by saying that "if you like excitement, conflict, 
t-H . and controversy, especially when nothing very serious is at stake, then you will love the history of quantization on 
curved spaces." ]l| The possibility of resolving the ambiguities experimentally for rigid body systems was the subject 
ON 1 of a conference discussion in 1957 among DeWitt, Wheeler, and Feynman The problem was also treated in Dirac's 
^ ■ famous text of 1967 S. Yet many years later we still seem far from attaining a consensus regarding quantization 
i[ on curved spaces. The "quasipermanent discussion among the specialists" ^] has by no means come to an end. In 
+Jj ■ fact, a recent paper Q laments that "in spite of all the successes of quantum mechanics and after years of efforts 
on quantization, we are still not absolutely sure about the correct quantization of as simple a system as the double 
pendulum." The subject also received early attention in the works of Cheng || and DeWitt |l. In this paper we 
do not attempt to review the extensive research on the subject of constrained quantization, and refer the interested 
reader to the literature ftm. Our approach to this problem is quite independent of all previous methods. 

The canonical quantization prescription has served us well in defining a straightforward and unambiguous way of 
quantizing a classical system in flat space. However, it immediately runs into problems where constrained dynamics is 
involved. It is easy enough to restrict the potential term to the constraint surface - the difficulty arises in treating the 
operator ordering ambiguities in quantizing the kinetic term K = — ^ -J~ (g 1 ^ ^fg ) . These ambiguities, at order 

h 2 , are proportional to the local Gaussian curvature of the constraint surface, which is the only available coordinate- 
invariant quantity. Thus for constrained systems with flat constraints (e.g. the simple pendulum or a particle confined 
to a one-dimensional curve), the ambiguity does not arise, whereas for systems with constant curvature (e.g. motion on 
a sphere) it leads only to a physically unobservable shift in the zero-point energy. In the latter case, the semiclassical 
Van Vleck approximation can be used to fix the "correct" energy offset j|] . In three dimensions this requires adding an 
R term to the Laplacian (R being the scalar curvature). For a two dimensional surface of constant negative curvature, 
the semiclassical approximation becomes exact (though only at distances large compared to the curvature radius) if 
the quantity added to the Laplacian is R/4. Even if one (somewhat arbitrarily?) requires that Van Vleck should be 
as close as possible to being exact for the constant curvature case, it is by no means clear that this prescription is 
valid in the generic case of varying curvature, where the ambiguities are most physically relevant. 

One might think that Feynman's path integral approach || would help to resolve these uncertainties. However, as 
remarked by Schulman in this context, "there is no free lunch." In the path integral method, the ambiguity arises in 
choosing how to incorporate the metric into the kernel and in deciding at what point in the infinitesimal time interval 
to evaluate functions of the metric [jlJQ. Two commonly used kernels give answers that differ by an effective R/6 
contribution to the Laplacian (and naturally enough both differ from the canonical quantization result). These issues 
are discussed at some length in Ali's recent paper ||, and here we will not go into the details. 
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Given this history, we have decided to attack this old problem using a somewhat different, and in our view, more 
physically motivated approach. The principle governing our view is that arriving at a unique quantization is essential 
only in real systems. That is, quantization should be physically unambiguous. For example, quantization of the usual 
rigidly constrained classical double pendulum is not a physical problem, because the requirement that the pendulum 
lengths remain exactly constant is not realizable. To make the problem as close as possible to physical, the constraints 
should be imposed through a limiting procedure, fluctuations remaining physically allowable along the way. We call 
this the "limit quantization," as opposed to direct quantizations which do not pass t hro ugh a limiting procedure. 
This procedure is carried out in Sec. |l| using a Hamiltonian formalism, and later in Sec. |lll| from a path integral point 
of view. The limit quantization embodies the idea of classical constraints arising as the limit of ever larger restoring 
forces normal to the constraint surface, an idea discussed in Arnol'd, for example H (see also Lanczos JnJ). In this 
way the constraint coordinates appear as part of the full dynamical system in the usual Euclidean metric. As the 
constraint forces are increased, the constrained degrees of freedom acquire a very high frequency, and their actions 
remain adiabatically constant as the slow variables move. This however does not imply that the fast variables (which 
we here call the r— variables) are energetically decoupled from the slow (9—) variables. We discuss the classical and 
semiclassical issues of constrained dynamics further in Sec. where from still another perspective we will see that 
the direct quantization is intrinsically ambiguous at order h : different limiting constraint potentials give the same 
classical constrained motion but different quantum dynamics at order h 2 . 

The uncertainty principle requires that in any real physical system there will be quantum fluctuations around the 
surface of constraint. By examining the effects of these fluctuations in the limit of ever tighter constraints, we find 
unavoidable ambiguities in the quantum dynamics along the surface of constraint. Our claim is that these ambiguities 
are an inherent part of any reasonable quantization procedure, and need not be more bothersome in principle than 
the Aharonov-Bohm effect, which shares the property of identical classical dynamics giving rise to different quantum 
dynamics. As the Aharonov-Bohm effect illustrates, there is no good reason why a quantum system should be uniquely 
determined by its classical (or even semiclassical) limit, and the richness of the new phenomena that can occur in 
physical quantum systems can well be viewed as a positive aspect of the quantum theory rather than as a shortcoming. 

In any case, we hope our contribution to the literature will be a valuable one, in spite of this concluding warning 
by Schulman: "Additional papers on curved space path integration are legion, and we do not attempt to list them all. 
Some are correct; some are less correct. Some have original features; some are less rich in this praiseworthy property." 

II. QUANTUM MECHANICS OF CONSTRAINED SYSTEMS 

A physically natural way of defining a constrained dynamical system in classical mechanics is through a limiting 
procedure where a strong attractive potential forces the system at any fixed energy to live closer and closer to the 
constraint surface |)| . The effective classical dynamics obtained in this limit turns out to depend only on the intrinsic 
properties of the constraint surface (e.g. the curvature) and not on the details of the constraining potential chosen 
on the embedding space. We will use an analogous procedure in quantum mechanics and consider the quantization 
of a constrained classical system as the limit of quantizations of classical systems with large potentials away from the 
surface of constraint. This approach is in accord with our intuition about physical constraints. It is also appealing 
because once a specific squeezing potential is chosen, the resulting quantum mechanics is entirely unambiguous, and 
can be obtained equivalently via a canonical or path integral method. Ambiguities will nonetheless arise in this 
approach, but are due to the freedom in selecting a constraining potential. The effective quantum dynamics in the 
squeezing limit will have corrections at order h 2 which depend on the details of the way in which this limit is taken. 
These ambiguities are analogous to operator-ordering effects in canonical quantization approaches, and to ambiguities 
which arise in the choice of kernel (particularly in the choice of evaluation points for functions of the metric) in path 
integral methods. However, the physical ambiguities we find are considerably more general in that they need not be 
functions only of the intrinsic properties of the constraint space. Even for a flat constraint space a large set of possible 
quantizations are available, each having its own physical realization in an embedding space, but all having the same 
classical and semiclassical limits. The existence of such quantum ambiguities in constrained systems all leading to 
the same classical physics in our view overshadows any attempt to resolve the operator ordering questions and terms 
involving the curvature. 

We begin with a Lagrangian 

L=\g ij q i q j -V{q)-V X {q) (1) 
on the full coordinate space, with flat metric gtj and potential V (we take the mass to be unity throughout). The 
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additional potential V\ enforces the constraint in the limit A — > oo. Classically we require V\ to have the same value 
everywhere on the constraint surface S, for each A. We now transform to the Hamiltonian 

H=^gVp iPj + V(q) + V x (q), (2) 

where pi = dL/dq l = g^q 3 is the momentum coordinate and g tJ is the inverse of gij. Canonical quantization now 
requires that we replace g^PiPj — p • p with the operator — ?i 2 V 2 , which in coordinate system q has the form 

where g is the determinant of the metric (i.e. the square of the volume element). Notice that there is no operator 
ordering ambiguity here because the full metric g^ is flat. 

Now locally, near some region of the constraint surface S we may separate the coordinates q into "slow coordinates" 
9 parametrizing S and "fast coordinates" r orthogonal to the surface (see Fig. (Q)). We will be looking for an effective 
theory of the coordinates 9 when the constraining potential (r) allows only small fluctuations in r (small compared 
to the length scale associated with the physical potential V, the curvature scale of the constraint surface S, and the 
scale on which V£ varies as a function of 9). The quantum wavefunction on the full space can be written as 

¥(g) = (0, r) = <h(P)*QR(r) + 01 ( )*i 8) (4) 



where (j>o, 4>i, ■ ■ ■ are arbitrary functions of the slow variables. "^grM is the ground state of the fast variables 
treating 9 as a fixed parameter, \tf^(r) is the first excited state, and so on. More precisely, *gr, ^i, etc. arc the 
eigenstates of the fast Hamiltonian 

B r = -~^V9^n^(T), (5) 

where i and j are summed over the fast variables only and the metric is evaluated at a fixed value of 9. Note that g 
here is the determinant of the full metric on the (r, 9) space. 

Now we want to consider the effective Hamiltonian H c g acting on the slow wavefunction 4>$(9) which multiplies 
the ground state of the r variables. We will argue towards the end of this analysis that this effective theory is in fact 
unitary in the limit A — > oo, as the transition probabilities connecting the ground states to the excited states of the 
fast variables disappear in the constraint limit (due to adiabaticity) . 

We have mentioned previously that in order to obtain the correct classical motion in the constraint limit, we must 
make V? (r = 0) be independent of 9 and have the potential increase away from r = 0. In order for the quantum 

mechanics to have a sensible semiclassical limit, we must further enforce that the energy of the ground state £qr, 
defined by 

^=4*1. (6) 

must be 9— independent, for any A, and the energy of the first excited state must be well separated from it, with 
the separation growing with A. This can be implemented in a number of ways, but two particularly simple scenarios 
present themselves. In the first we take V\(q) = Au(g), where v is a smooth function vanishing on the constraint 
surface and having its minimum there (but with a non-singular second derivative matrix with respect to the fast 
variables). Near the constraint surface, the constraining potential will have the harmonic oscillator form 

Ve(r)=r i A«? r i + 0(r 3 ), (7) 

where the matrix function is smooth. The ground states will to first approximation be harmonic oscillator 

ground states with spatial extent of order ft 1 / 2 /A 1 / 4 and energy -Egr of order TiX 1 / 2 . As discussed above, the coefficient 
of this energy must be ^-independent, to prevent infinite effective forces in 9 from arising in the A — > oo limit. The 
matrices A^ s > can easily be adjusted to satisfy this condition. In the case of more than one constraint variable, this still 
allows for much freedom in the function A^ and the resulting eigenstates (energy can flow from one fast degree 
of freedom to another as a function of 9, as long as the total energy remains constant). One may choose to adjust 
(some of) the anharmonic parts of the potential as well but this is not really necessary or natural. The quartic term 
will lead to order fi 2 (A-independent) corrections to the ground state energy and thus (if these anharmonic corrections 
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vary with 9 which in general they will) to the effective potential in 9. [] We will later see that this effect is of the 
same order as other terms that are encountered in the effective potential. We might also mention at this point that if 
in finding the ground state of r at a given 9 we used normal coordinates centered at (6, r — 0), gij = 6y + 0(r 2 ) for 
fixed 9, and neglected the 0(r 2 ) curvature correction, we would also make an error of order h 2 in the ground state 
energies and therefore in the effective potential for 9. 

A second scenario involves a hard-wall potential with the cross-section of the allowed region shrinking towards S 
as A — > oo. Thus we may instead take 

V°{r) = 0, aW(Ar) < 1 (8) 
co, aW(Ar) > 1, (9) 

where the function aS 6 ' defines a 9— dependent region shrinking as 1/A in which the fast eigenstates live. The energy 
of the ground state is of order U 2 \ 2 and can again be adjusted to be 9— independent while allowing for substantial 
freedom in the shapes of the ground states. In the following we will focus on the original (smooth potential) scenario, 
but this one could be used to obtain similar results. 

We now see what happens when we compute the effective Hamiltonian for the slow variables H c g(9, d/d9), defined 
by 

(^Iffcfflfo) = <*'|ff|*> (10) 

with 

*(0,r) = 0o(0)*g£(r). (11) 

The potential V{q) clearly induces an effective potential V(6,r = 0), up to corrections which disappear in the 

constraint limit, since V(9,r) varies slowly in r over the extent of ^Q^(r). The constraint potential V\, combined 
with the fast part of the kinetic term, gives, by assumption, a ^-independent effective potential which we can therefore 
drop. We are left with the slow (0) part of the full q— space Laplacian, as well as mixed terms involving one derivative 
each with respect to the slow and fast variables. 

In the Born-Oppenheimer spirit, consider first the action of the "slow" part of the Laplacian on the full wavefunction. 
We are interested in matrix elements of the form 

(^* GR | - Y^W {9ij ^W Mo * GR) ■ (12) 

The two derivatives can act on the slow wavefunction <f>o, the ground state \I/gr> or on the metric and its determinant. 
Consider terms where both derivatives act on the ground state or the metric. After integrating over the fast variables, 
this leads to (cf/ \ — ?i 2 /(#)|0o), an effective potential scaling as h 2 . [] (Note that both the metric and the ground 
state are changing over typical scales of order unity in 9.) Furthermore, this effective potential is A— independent for 
large A, since "Jgr maintains its functional form (that of a gaussian plus higher order corrections) and simply shrinks 
towards r = as A — > co. 

Now suppose one of the derivatives acts on the metric or on ^gRj while the other acts on the slow wavefunction. 
We then obtain 

(^|- fi2ftl (^+^-IW. (13) 

where h % are real functions of 9 (because the ground state ^gr can be chosen to be real, and the metric is real as well) . 
Taking the hermitian part we again obtain an effective potential of order h 2 . Finally, if both derivatives act on the 
slow wavefunction </> , we obtain the usual kinetic term for the slow variables, including order Ti 2 terms proportional 
to the curvature of the constraining surface, as discussed in the preceding section. This concludes the discussion of 
the pure 9— derivative part of the kinetic energy, and we move on to the mixed terms next. 



1 For the purpose of dimension counting, every power of r in an expression contributes a factor of h 1 ^ 2 /A 1 / 4 to the expectation 
value whereas d/dr acting on *1/gr (but not the metric) contributes A 1//4 /H 1 ^ 2 . Derivatives acting on the "slow" wavefunction 
4>o{0) and the metric g, as well as d/d9 acting on ^Pgr all produce factors of order unity. 

2 The derivatives acting on the metric only will give the usual curvature-dependent corrections to the effective potential. Terms 
coming from derivatives acting on the ground state are of the same order and will give rise to additional corrections. 
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If the d/dr derivative acts on the metric, we obtain a contribution of order ft 2 to the effective potential as before 
(all derivatives are of order unity). If d/dr acts on \I/gr while d/d9 acts on g or the wavefunction <fio, we obtain 

h 2 (* GR \b e (r)-?-+h.c.\* GR ), (14) 

where b is a slowly varying real function. So to this order in ft only the commutator of b with d/dr survives, and once 
again the resulting contribution to the effective potential is of order ft 2 . The remaining term is the one where both 
derivatives act on the ground state; it has the form 

d d 

h 2 (y GR \c(e,r)— — \y GR ). (15) 

Since 4'gr('*) is even in r, the slowly varying function c(9,r) must be expanded to first order in r to get a non- 
vanishing contribution. We then obtain yet another contribution to the effective potential of order ft 2 (since r 
acting on ^grM is of order unity). 

Finally, we argue by adiabaticity that there is no mixing between 4>o(9) and the excited wavefunctions (j)i(9) 1 etc. 
This follows because the splitting between ^ GR and the excited states of r scales as ft\/\ in the constraint limit (or 
as ft 2 X 2 in the case of a hard wall potential), whereas the matrix elements do not grow with A. So assuming a finite 
amount of energy in the initial (/)q(9) wavefunction, transition probabilities vanish in the A — > oo limit, and we obtain 
a probability-conserving effective quantum mechanics. The effective Hamiltonian is 

= ~^W^^ + V{6) + h2VM • (16) 

This is the key result of this paper. V c s is essentially an arbitrary smooth function, depending in a complicated 
way on both the curvature properties of S and on the constraining potential around S (together these determine the 
ground state wavefunctions). ^fg can be taken to be the square root of the determinant of either the full metric or of 
the metric restricted to the constraint surface. Any difference between these expressions can be absorbed into V e s. 



III. THE LAGRANGIAN FORMALISM 

A path integral analysis also shows that quantization on the constrained surface is ambiguous at order ft 2 . We 
begin with the propagator on the full coordinate space. Because the space is flat, this is given by the usual Feynman 
prescription J 2?(q(r)) exp(iS(q(T))/ft) where S(q(r)) = J L(q(r),q(T),T)dT. The Lagrangian L on the full space is 
given by Eq. |l| in Section Q we follow the notation introduced there. We extract the effective dynamics for the 
"slow" variables 9 by performing a trace over the "fast" variables r in the adiabatic limit. Specifically, we calculate 
(9 f\p(t)\9 ) where pit) = Tr r (W(t)) is the trace over the fast variables of the density operator on the full space. The 
time t is taken to be small on the time scales of 0-motion but r may undergo many oscillations during this time. 
Noting that the evolution of the full density operator is given by W(t) — e ~ lHt / h W(0)e tHt / h , and inserting identities 
in the form of complete sets of states, we have 

(9 f \p(t)\9 ) = f dRdr'dr''d9 / d9 / \9 f R\e~* m/h \9'r')(9'r'\W{Q)\9"r")(9"r"\e im/n \9 R) . (17) 



We take the initial density operator to be the "product" state I^qr.) ® l^)> wriere the initial state of the fast variables 
is the (purely real) ground state ^^(r), with 9 regarded as a fixed parameter (see Section ||). 4>{9) is an arbitrary 
initial wavefunction of the slow variables 9. Replacing the two propagators by path integrals, we then have 

(9 f \p(t)\9 ) = [ d9'd9" [ m 9, V9(t) [ m 6 ° V9{r)cp{9'W{9")T (9{t),9{t)) , (18) 
J Je(o)=e' Je(o)=9" v 7 

where 

■F(0(t),0(t)) = J dRdr'dr" 

r r(t)=R r r(t)=R 

x / Vr{r) / Vr(Ty s(1) / h ~ tS(2) / h ¥° R \r")¥°2(r'). 

Jr(0)=r' Jf(0)=r" 

(19) 
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SW — S(r(r); #(r)) and = S(f(T);9(T)) are the actions along the respective paths. In performing the path 
integrals over the fast variables, 9(t) (9(t)) is to be treated as an external forcing function with the property 9(0) = 
8', 9{t) = Of (resp. 0(0) = 9", 9{t) = 9 a ). 

We evaluate T using stationary phase with respect to Ti on the path integrals, trace, and the integrals over the 
intermediate fast variables. In fact the stationary phase approximation becomes exact in the adiabatic limit: we shall 
now show that the errors due to stationary phase evaluation are of order ft^/A 1 / 4 , so they vanish as A — > oo. Recall 
that 



J g{x)e if ^ h dx « sp 



I 2nih 
f"{x n ) 



g{x r , 



i + o(h 



1/2 I'"{ x n) 



(20) 



where the sum is over stationary phase points x n satisfying f'(x n ) = 0. Stationary phase evaluation of path integrals, 
although different in the details, scales in the same way. f(x) corresponds to the action S. The stationary paths are 
the classical paths and we now argue that S c \ ~ 0(A 1//2 ) to leading order in A. The leading behavior of the action in 
the fast variables arises from the kinetic term in r and the constraint potential: L r ~ gy-f i r- 7 /2 — V\(r), where the 
sum is over fast variables only. This Lagrangian corresponds to the fast Hamiltonian of Eq. |^ discussed in Section [n|, 
which gives harmonic motion for r at least throughout the range of the fast variable ground state. Because the initial 
fast state is this ground state, the leading behavior of e iS ' n gives 



dr' / Vr(r)e 



iS 



(I,/ft 4 e aV) 



(21) 



where we note that -Egr ~ 0{hy/X) is required to be 0-independent to avoid infinite torques acting on the slow 
variables (Section ||). This shows the leading behavior of the action on the ground state is simply to multiply it 
by an evolving phase. The fast variable r(t) does not stray from the domain of faithful harmonic approximation to 
the constraint potential. Subleading terms are a factor of 1/y/X smaller. Thus the dominant behavior of the action 
and its derivatives is simply that of a harmonic potential of frequency O(vA). The errors in the stationary phase 
evaluation of the path integrals are of order h 1 ^ 2 /X 1 ^ 4 , vanishing in the adiabatic limit. The errors in using stationary 
phase to compute the trace and integrals over the intermediate variables r' and r" scale in the same way. In the 
harmonic approximation to the action and the ground state it is readily seen that the stationary phase points are 
= r" = 0, the constrained value of the fast variable. Corrections to this approximation due to subleading 



R = r' = 
terms vanish in the adiabatic limit, 
We have then, 



A 



T(9(t),9(t)) 



d 2 In* 



dr 1 '- 



2e(i) 



d 2 S. 



dRdr 1 



1/2 



-1/2 



GR 



d 2 s 



(2) 



dRdr 



' d 2 ln^R 



-1/2 



dr" 



1/2 



(22) 



where the derivatives are evaluated at R = r' = r" = 0. S^j (resp. Sf?') is the action along the classical path starting 
and ending at the constraint value r = subject to the Lagrangian of Eq. with 9(t) (resp. 9(t)) treated as an 
undetermined forcing function. In the adiabatic limit, we shall now show that the action exponent is just the action 
for the reduced slow variable system on the constraint surface, i.e. the action we would have written down had we 
begun in the reduced space. The kinetic term in r together with the constraint potential give a ^-independent term in 
the exponent as discussed above. This, being a constant energy shift as far as 9 is concerned, can be neglected. The 
mixed kinetic term goes to zero: we may expand gij(r, 9) about r = and take functions of 9 out of the time-integral 
as they are slowly varying functions of time: 



9lJ (r,9)r l ¥ dt = g v (0,9)& f l dt 

Jo 

+ 0(n 3 / 2 /A 1/4 ) , 



drP J 



(23) 



where the derivatives of the metric are evaluated at r = 0. The first term on the right-hand side is zero, since the 
integral gives r l (t) — r l (0) = due to the stationary phase condition on the endpoints. For p — i, the integral in the 
second term is zero for the same reason (the integrand is d(r l2 )/dt). For p ^ i the integral averages to zero, because 
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different directions of the fast variables generically oscillate at different frequencies. The higher order terms in the 
expansion vanish in the adiabatic limit (they involve at least two powers of r but only one power of r so scale at least 
as (%/X 1 / 2 )(h 1 ' 2 X 1 ^)). We are thus left with the kinetic term in the slow variables only, with the metric evaluated 
at the constraint value for r. 

Inserting J-(8(t), 8(t)) into Eq. |l8| gives the reduced density matrix at time t. This factorizes into a part involving 
(9,8 f, and 8(t)) and a part involving (8 1 ,8 , and 6(t)). This implies that the reduced density matrix factorizes: 
(8f\p(t)\8 ) = (9f\8(t))(0(T)\6 o ), i.e. the final 0-state is pure and so we may describe it in terms of a wavefunction. 
Adiabaticity has thus uncoupled the fast and slow degrees of freedom. We have 

cf>(9 f ,t) = [dd [ V8(r)A(8(r))e lS l h 0(8,0), (24) 



where 



A(8(r)) 



dr 



GR I 
; r=Q 



dRdr 



1/2 

(25) 



S = is the action function for the 8 variable on the constrained surface: S = J Ldt where the Lagrangian is 
L = ^gijd 8^ — V(8), with gij now the metric on the curved space defined by the constraint surface and the sum is 
over the slow ((9-)variables only. 

We may compare this to the expressions discussed in Ali, where one works in the reduced space from the start. 
Different Feynman kernels are postulated to attempt to account for the curvature of the constrained surface |l],f|,[| ■ 
There, for infinitesimal time t, 

<f>(9 f ,t) = I dn(9)G(6,9 f )e iS / n (f>(9,0). (26) 



dQ(9) = \J g(8)dd is the volume element at 9. Candidates for G(8,8f) which are often considered in the literature 
include the identity operator and g(8)~ 1 / 4 D 1 / 2 g(8 j) 1 / 4 , where g is the determinant of the metric on the curved space, 
and D is the Van-Vleck determinant (D — det(—d 2 S/d8d8f)). The different choices give rise to Hamiltonians which 
differ at order ft 2 by a certain fraction of Ti 2 R. Ali demonstrates this by choosing locally normal coordinates £ to 
evaluate the integrals and considering infinitesimal time so that the exponent and prefactor can be expanded in £. The 
action exponent is then at lowest order (in £) quadratic, and the resulting Gaussian integrals are readily performed. 
Only the even order terms in the expansion of the prefactor about the initial value 8 contribute. The t — > limit is 
taken and an effective Hamiltonian can be extracted from the resulting differential equation. The expansion of the 
potential exponent and the initial wavefunction give the Hamiltonian — 7i 2 V 2 /2m + V, and it is the quadratic term 
in the expansion of \/ g(£)G(£) that gives rise to the Ti 2 R discrepancies in the effective Hamiltonian. (Higher order 
terms give corrections at higher order in h.) 

The key point is that we may apply the same manipulations to our Eq. |^, where A(8(t)) plays the role of the 
prefactor y/ g(9)G(8). Transforming to locally normal coordinates in which there are no linear terms in the action 
exponent, and expanding A(9(t)) about 9 leads to finite Ti 2 corrections in the effective Hamiltonian as in Ali's 
approach described above. However, it is important to recognize that our corrections emerge from reducing the full, 
unambiguous coordinate space evolution down to the reduced curved space by taking the constraint limit, whereas 
Ali's corrections arise from using a modified Feynman kernel for curved space. We observe that our corrections depend 
not only on the curvature but also on the constraint potential. These conclusions are in accord with the results of the 
previous section and provide additional insight into the ambiguity and connection with the literature. 



IV. CLASSICAL AND SEMICLASSICAL PERSPECTIVE 



Finally we provide a third, classical and semiclassical perspective on the ambiguities inherent in quantizing con- 
strained systems. The classical point of view on constraints as limits of strong restoring forces is instructive, and nicely 
complements the quantum discussion. In classical mechanics, adiabaticity plays an important role. Adiabaticity of 
classical actions also played an important role in the old quantum theory. Here we will see that the semiclassical 
approach confirms the inherent order h 2 problem of quantizing a constrained system. Again the plan is to start 
with the system in the full Cartesian space, transform to curvilinear coordinates and the Laplace-Beltrami kinetic 
energy operator with no ambiguity, and begin imposing stiff force constants about fixed values of one or more of the 
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coordinates. The result in the limit of infinite restoring forces is constrained dynamics on a Riemannian manifold. 
Again an ambiguity will arise because we can arrive at identical classical constrained dynamics in ways which differ 
quantum mechanically. 

The idea that classical constraints are the limit of merely stiff degrees of freedom is a natural one, though not often 
discussed. Arnol'd includes a brief treatment of this in his famous book in a chapter devoted to Lagrangian dynamics 
on manifolds |§. Arnol'd states that in the limit A — > oo the classical dynamics satisfies Lagrange's equations of 



motion 



d ( dL* \ _ f dL. t 



(27) 



where 

L*{{0}, {8}) = 21^=0^=0 - V\ ri=0 , (28) 

where T is the kinetic energy and V any non-constraint potential energy of the full system. However there is an 
important difference for the present purposes: we need to consider very large energies stored in the constrained 
coordinates, corresponding to the quantum zero point energy climbing to arbitrarily high values as A — > oo. This we 
can do, however, using the theorem on the adiabaticity of classical actions under slow variation of a parameter. 

The issues which confront the classical approach to constraints in the high energy constraint limit are 1) energy 
stored in the fast variables generally varies as a function of the slow variables, even though actions are constant, 
thus raising the spectre of large energy exchanges with the slow degrees of freedom, 2) the possibility that the slow 
coordinates are chaotic, and 3) the possibility that the constraint coordinates can exchange energy among themselves, 
perhaps chaotically. 



A. One dimension 



In the simple case of one dimensional motion with an adiabatic change of one or more parameters, the variation 
of the action decreases exponentially as the rate of variation of the parameters decreases [[ll] . However as the phase 
space diagram in Fig. (|[a) shows, the energy generally changes significantly as the parameters change, as measured 
by the contemporary value of the Hamiltonian. For example, in the case of a one dimensional harmonic oscillator, 
suppose we slowly adjust the frequency u>(t) = A(t)wo toward higher values, beginning with some energy Eq in the 
oscillator: 

H=~i? + ~\{t) U tf. (29) 

Since the action is constant, the area of the ellipse enclosed by the trajectory is constant, but this means that at q = 
the momentum must increase. However at q = 0, H = p 2 /2, so the energy increases as u(t) increases adiabatically. 
Specifically, with uj oc vA, we have p oc A 1 / 4 and E cx A 1 / 2 . 

Now suppose we have only one constraint coordinate r. The situation is depicted in Fig. The slow 9— coordinates 
act as adiabatic parameters. Whatever their motion, in the limit of large, fixed A the action of the r— motion is 
preserved during the course of the 9 motion. Clearly if we are to avoid any of the large zero point energy of the 
r— coordinates finding its way into the 9— coordinates we must restrict the confining potential V(r; 0±, 9 2 , . . .) to vary 
with 9±, 62, ■ ■ ■ in a way that maintains fixed energy for fixed action. This condition does not fix the potential uniquely, 
but is possible for a wide class of potentials, for example satisfying 

E X {I) = A 1//2 (ciJ + c 2 I 2 + ■■■) = A 1 / 2 e(7) (30) 

with 



e(ft/2) = eo (31) 
independent of A. The action / is defined as usual as 

I=^j>p r {E,\)dr. (32) 

An example showing two potentials (or rather their phase space contours at identical actions and energies) satisfying 
this requirement is shown in Fig. (||,b). The area enclosed by the dashed and solid contours is the same and the 
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momentum at q — (where we assume V = for both potentials) is also identical, showing the energy is the same 
for the two potentials at that action. Any potential satisfying Eqs. ( fjol) and ([H]) at fixed values of the 9— coordinates 
is an allowed constraint potential at that point in 9 space; the independence of Eqs. (||(]) and ( |3l| ) under changes in 
the 9— coordinates must then be arranged to ensure energy remaining fixed in the r— variable. 

The semiclassical estimate of the energy in the r— coordinate is thus independent of the 9— coordinates, since that 
energy is just given by the classical energy at action / = Ti/2; this is what we set constant by definition. However it 
is not the semiclassical energy we want, but rather the exact quantum energy. Now the key fact is that the exact and 
semiclassical energy differ in order Ti 2 . This follows from the derivation of the WKB energy equation, which is an 
expansion in fi. Thus, the freedom to choose among potentials which lead to identical classical constrained motion as 
A — > oo leads unavoidably to a quantum ambiguity at order h 2 , in agreement with the discussion in the previous two 
sections. This is another confirmation of the view forwarded there, and is the main result of this section. 



B. Chaos in {8} Variables 



The issue of chaos in the {9} variables is easily disposed of: it makes no difference, since in the limit A — > oo the 
frequency of the r— variable is far from reach of the slow 9— coordinates, which just act as slow parameter changes on 
the r— variable motion. The fact that the changes of the 9— coordinates are not quasiperiodic in the chaotic case is of 
no consequence to the adiabaticity of the r— coordinates. 



C. Resonance or Chaos in {r} Variables 



Interaction between several r— variables can be avoided by independent control of as many A— parameters as 
r— coordinates. With this control we can force the frequency ratios uii/wi+i — > even as each u>i — > oo. By this 
ruse we avoid any low order resonance leading to energy transfer between the r— variables. The concern that high 
order resonances must always exist is softened by the fact that even if the external control parameters (9— variables) 
are effectively fixed and the high order resonance is allowed to act, the energy exchange (resonance width) is expected 
to be exponentially small in the winding number (frequency ratio) for smooth potentials, as it depends on Fourier 
coefficients of the order of the winding number. In effect, good actions are maintained in the r— variables as A — ► oo. 

In the opposite extreme of full chaos in the r— variables, we appeal to Liouville's theorem applied to the invariance 
of the total phase space volume enclosed by the energy surface in the r— motion in the limit that the 9— variables 
are slow enough to be considered to be adiabatically separate external parameters. In addition, we need Hertz's 
theorem, which states that in an ergodic system the energy shell maps onto another energy shell in the adiabatic 



limit of changing parameters 12 15J. We may write E — E^(V), where V is the phase space volume. As long as 
the constraint potentials are tuned to preserve this energy as a function of the 9~ variables, i.e. E^gj(V) = E^giy(V) 
for all {9'} and fixed V, there will again be no forces on them from the constraints. Having arranged the potentials 
in this way, there remains the question of the semiclassical energy for the chaotic dynamics. This can be given by 
periodic orbit theory fl4|| . Somewhat paradoxically, periodic orbit theory has been shown to give excellent results for 
low energy eigenvalues in a number of systems, breaking down at higher energies. We need the lowest (zero point) 
energy, which at least formally has an error of order h 2 , so once again we reach the same conclusion about the order 
h ambiguity in the quantization procedure. 



V. CONCLUSION 



We have shown from three different points of view that the ambiguities which arise in the quantization of a 
constrained system are physical in nature and cannot be resolved by finding the "right" mathematical formalism. 
Considering the constrained surface as a reduced sub-space of a larger system, we have found that terms of order h 2 
inevitably arise in the effective description. These terms depend on the details of the forcing potential as much as on 
the intrinsic geometry of the constrained surface, but there is no ambiguity in the classical or even the semiclassical 
dynamics. Thus we obtain a large set of equally valid "quantizations" of the same classical system, each corresponding 
to a particular constraining process. In any physical situation, one such constraining process is the correct one, and 
the ambiguities disappear. In the absence of detailed knowledge about the nature of the constraints, additional 
"quantization conditions" must be selected before one can speak of quantizing any classical dynamics on a curved 
space. 







Our discussion has been cast in terms of holonomic constraints (which can be written in terms of the vanishing of 
one or more functions of the coordinates) but the idea of reaching rigid constraints via a limiting process applies also 
to nonholonomic constraints, which can be couched only in terms of nonintegrable differentials. 
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Figure Captions 

FIG. 1. Schematic showing the fast motion along r and slow motion along 61, #2 coordinates for finite confining forces (finite 
A). 

FIG. 2. a. An area preserving (adiabatic) change of Hamiltonian parameters generally changes the energy, as indicated by 
the increase in \p\ at q = 0, where the potential V = 0. b) Some slow deformations of the Hamiltonian do leave the energy 
unchanged. 
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